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Surface codes have emerged as promising candidates for quantum information processing. Building 
on the previous idea to realize the physical qubits of such systems in terms of Majorana bound states 
supported by topological semiconductor nanowires, we show that the basic code operations, namely 
projective stabilizer measurements and qubit manipulations, can be implemented by conventional 
tunnel conductance probes and charge pumping via single-electron transistors, respectively. The 
simplicity of the access scheme suggests that a functional code might be in close experimental reach. 

PACS numbers: 03.67.Pp, 03.67.Ac, 03.67.Lx, 74.78.Na 


Introduction .—In recent years, surface codes have es¬ 
tablished themselves as a potent platform for universal 
quantum information processing. The basic idea of a sur¬ 
face code (when operated as a so-called stabilizer code) 
is to implement few ‘logical’ qubits - the actual carriers 
of quantum information - via the correlation of a large 
number of physical qubits What at first sight ap¬ 

pears to be a redundant scheme offers a number of pow¬ 
erful advantages: (i) a degree of tolerance to errors or¬ 
ders of magnitude higher than that of other approaches, 
(ii) the possibility to implement a code as a compara¬ 
tively simple two-dimensional (2D) layout of cells cou¬ 
pled by nearest-neighbor interactions, and (iii) the fact 
that essential code operations, including error tracking 
without need of active error correction in Clifford oper¬ 
ations and/or memory access, are controlled by classical 
software. (While non-local stabilizer codes may allow 
for even higher error thresholds [6], we here focus on lo¬ 
cal surface codes.) These are highly attractive features 
which, accordingly, come at a hefty price tag: a large 
number of physical qubits is required even for modest 
operations. In particular, logical non-Clifford gates (e.g., 
the T gate) are required for universality, whose fault- 
tolerant implementation would require magic state dis¬ 
tillation ElIH]. Under these conditions. Ref. [5] estimates 
that about 10^“^ physical qubits are needed to encode 
a reasonably fault-tolerant information qubit, implying 
that about O(IO^) physical qubits are needed to run, say, 
serious factorization algorithms for integers with 0(10^) 
decimals. Achieving maximal simplicity in the implemen¬ 
tation and in the access of individual qubits will there¬ 
fore be a decisive factor in advancing from 0(1) qubits 
to functional systems. In this Letter we argue that semi¬ 
conductor Majorana hardware layouts offer striking and 
so far unnoticed advantages in this regard. 

At this point, two major platforms for the realization 
of surface codes are under discussion, Josephson junction 
arrays, and Majorana bound state (MBS) networks, re¬ 


spectively. The Josephson junction architecture builds 
on physical qubits that are an experimental reality, and 
impressive progress towards the generalization to qubit 
assemblies has been made recently HHn]. By contrast, 
not even the building blocks of an MBS qubit have been 
implemented so far. However, few years after the pre¬ 
diction [IMS] of semiconductor MBSs, and the obser¬ 
vation of first experimental signatures m, a now avail¬ 
able second generation of topological nanowires features 
robust proximity coupling to an adjacent superconduc¬ 
tor m and sharply defined proximity gaps [IHlIini. It 
stands to reason that this sets the stage not only to the 
clean isolation of midgap MBSs, but also to their cou¬ 
pling by tunneling bridges m, which will be instrumen¬ 
tal to the definition of physical qubits based on MBS ring 
exchange m- 

A principal advantage of the semiconductor platform 
is that it can be accessed in terms of single-step protocols 
which do not require ancilla qubits [22]. (This is to be 
compared to a five-step projective measurement requiring 
microwave resonators and one ancilla qubit per physical 
qubit in approaches based on bosonic implementations of 
quantum storage laiMii) However, the concrete Ma¬ 
jorana code readout procedures suggested so far rely on a 
combination of gate electrode operations, microwave irra¬ 
diation, and subsequent measurement of the phase shifts 
in the transmitted photon field in a topological insulator 
platform [22], or on magnetic interferometry on the single 
qubit level within a nanowire setting m- Both schemes 
involve extensive hardware overheads which raises ques¬ 
tions regarding scalability. 

In this Letter, we show that MBS codes can be accessed 
in much simpler terms via ordinary tunnel electrodes and 
similarly established components of semiconductor quan¬ 
tum electronics that can be integrated into the 2D code 
structure itself. Within this framework, the system is op¬ 
erated by tunnel conductance measurements and single¬ 
electron transistor (SET) manipulations in combination 
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Figure 1. Structure of a surface code built on MBSs [red dot in 
panel (a)], few of which define a Majorana-Cooper box [(b)]. 
Ring exchange of MBSs via tunneling bridges [(c)] defines a 
physical qubit through the eigenvalues ±1 of the plaquette 
operator On, i-e., the product of all eight Majorana operators 
in the loop. An assembly of qubits [(d)] can encode informa¬ 
tion qubits, e.g., the unit indicated by a dark shaded area. 
For further discussion, see text. 

with classical software protocols. Importantly, no radia¬ 
tion fields nor probing magnetic fluxes are required. 

Semiconductor-based Majorana surface code .— To 
start with, let us consider the schematic blueprint of a 
Majorana code depicted in Fig. Its elementary build¬ 
ing blocks are MBSs forming at the terminal points of 
spin-orbit coupled (e.g., InAs) nanowires in a magnetic 
Zeeman field and proximitized by an s-wave supercon¬ 
ductor [HHnmsHii], cf. the red dot in Fig. [^a), where 
the wire and the superconductor are indicated by dark 
and light gray, respectively. We consider sufficiently long 
wires such that the MBSs, which are described by anti¬ 
commuting operators 7 ;l ,2 = 7 i ,2 ~ represent 

zero-energy states. Pairs of wires (a = 1, 2) are contacted 
to the same floating mesoscopic superconducting island, 
see Fig. Bb), such that the entire ‘Majorana-Cooper box’ 
[231127] has a finite capacitance. Each box then has a 
charging Hamiltonian He = Ec{M — ng)‘^^ where M is 
the total number of electrons of the island relative to a 
backgate parameter Ug , and Ec the single-electron charg¬ 
ing energy. Assuming that Ec defines the largest energy 
scale in the problem and that Ug is close to an integer 
value. He will enforce charge quantization on the box. 
The resulting parity constraint, na=l, 2 ^ 7 i 72 = ±1, re¬ 
duces the four-fold ground-state degeneracy associated 
with four MBSs to a two-fold degeneracy. This effective 
‘spin-1/2’ degree of freedom is nonlocally encoded by the 
MBSs and implies a ‘topological Kondo effect’ [21] when a 
single box is contacted by normal leads. Instead, we here 
consider a 2D array of Majorana-Cooper boxes whose 
MBSs are tunnel-coupled to semiconductor wires con¬ 
necting neighboring boxes I and see Fig. [^c). Assum¬ 
ing a near instantaneous transmission through the (short) 


tunneling bridges, we model the latter by weak static am¬ 
plitudes tiif with \tiij Ec- The corresponding tunnel¬ 
ing Hamiltonian reads Ht = —+ h.c., 
where cpi is the fluctuating superconducting phase of box 
I [24H26]. 

Under these assumptions, the low-energy excitations 
of the 2D array correspond to minimal-loop structures 
involving the product On = 11^=1 ^f the eight 

MBSs, surrounding loop no. n of the network. The 
plaquette operators On are the so-called stabilizers of 
the system. Their defining feature is that they mutu¬ 
ally commute (two stabilizers share either zero or two 
MBSs with each other) and possess only two eigenvalues, 
spec (On) = ±1. This makes the stabilizers a primitive 
system of qubits, the backbone of the surface code. The 
corresponding effective low energy Hamiltonian 

Eq ~ ^ YOE^ 1J. ’ (^) 

n ^ 

contains the products of the four tunnel amplitudes con¬ 
necting the islands In along the nth loop, where the nu¬ 
merical prefactor follows from a Schrieffer-Wolff trans¬ 
formation [28], see also Refs. [2M3T] . The projection to 
the lowest charge sector in each box effectively removes 
any dependence on the phases Lpi in Eq. §■ Notably, the 
code works even though the coefficients Re(cn) in Eq. 0 
are uncorrelated random energies in general [32]. 

Operation principles. —Let us recapitulate the essential 
steps by which a local surface code is operated [IHS]. At 
the heart of the procedure is a continued sequential mea¬ 
surement of all stabilizers On- This operation serves to 
project the code onto a well-defined simultaneous eigen¬ 
state of the stabilizer system. As far as memory access 
or Clifford operations are concerned, occasional flips of 
individual stabilizers due to thermal or other fluctuations 
need not be actively corrected; they are simply recorded 
by classical control. This simplification follows from the 
Gottesman-Knill theorem [34], and the reduced need for 
active error correction is one of the principal advantages 
of the scheme. Information qubits are imprinted into the 
code matrix either by ceasing the measurement at indi¬ 
vidual (or groups of) physical qubits, or by creating phys¬ 
ical defects through the removal of tunneling links [cf. the 
dark area in Eig.j^d)]. In either case, a ‘hole’ is punched 
into the system, and the binary eigenstates of a suitably 
constructed Z 2 Wilson loop operator mm surrounding 
the hole then define an information qubit. Einally, the 
states of both physical and information qubits must be 
adjustable via the controlled flipping of individual stabi¬ 
lizer states. With these operations in place, all quantum 
gates required for universal quantum computation can be 
implemented. Using magic state distillation [8], this can 
be done in a fault-tolerant manner [6]. Rather high mi¬ 
croscopic error rates are tolerable if only the ratio in the 
number of physical to information qubits is sufficiently 
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Figure 2. Surface code setup. Pairs of normal leads (indicated 
by vertical lines for clarity) are tunnel coupled to MBSs 71 and 
72 located on different boxes. The two-terminal conductance 
measurement provides information about plaquettes Oa and 
Ob, cf. Eq. Q. Using SETs in the external circuit, plaquettes 
can be flipped as discussed in the text. 

large. It is evident, then, that the efficiency of quantum 
information processing crucially depends on the availabil¬ 
ity of simple stabilizer readout and manipulation proto¬ 
cols, satisfying the criteria given in the Introduction. In 
the following, we discuss how such protocols can be im¬ 
plemented. 

Measurement layout .—We consider the setup in Fig.|^ 
under the assumption that normal electrodes are con¬ 
nected to individual MBSs. (For a more concrete dis¬ 
cussion of a device architecture, see below.) A tunnel 
coupling of amplitude Xj between the Major ana fermion 
7 j and a lead electrode (with fermion operator at the 
contact) is modeled by the Hamiltonian [26] . 

+ h.c., ( 2 ) 

where lowers charge by one unit on the box host¬ 

ing 7 j. Since jj anticommutes with the two plaquette 
operators On containing this MBS (and commutes with 
all others), flips precisely those two plaquette eigen¬ 
values but leaves all others untouched. It is then straight¬ 
forward to compute the current flowing between two tun¬ 
nel electrodes {j = 1 , 2 ) biased by a transport voltage. 
Such a current measurement will induce plaquette flips 
in general, with one important exception: if the leads 
are attached to neighboring MBSs located on different 
boxes, see Fig. both 71 and 72 belong to the same 
pair of plaquettes, On=A,B- The corresponding plaque¬ 
tte eigenvalues are thus flipped twice (and hence remain 
invariant) whenever an electron is transported through 
the system. 

To make these statements quantitative, we project out 
all higher box charge states by Schrieffer-Wolff transfer- 
mation, He + Ht + Y,j=i ,2 ^ where 

iJeff = a (^ +40 a+cbC>b) + h.c., (3) 

a = —32AiA2/(5t^2^C')5 and the presence of the plaque¬ 
tte amplitudes Cn=A,B of Eq. 0 indicates charge trans¬ 
fer processes around either loop A or B. The tunneling 


amplitude ^ = [b\ti 2 \‘^/{16Ec)]r] describes direct charge 
transport along the link 1^2 and exceeds the loop am¬ 
plitudes parametrically in {Ecltiy)'^. However, it is also 
proportional to a ‘detuning parameter’ r] [33], which for 
small deviation of the gate parameters Ug^i off in¬ 

teger values scales as 77 ^ Ang^iAng^ 2 - Specifically, for 
Aug^i = 0 , destructive interference removes the direct 
amplitude ^ ^ completely. Perturbation theory in i^eff 
then yields the tunnel conductance as 

= 47r^|(apz/iz/2 {qq + QaOa + QbOb + QabOaOb) , 

e^ / a 

( 4 ) 

where z/ 1^2 are the density of states of the leads, and Oa,b 
now represent the eigenvalues (± 1 ) of the respective pla¬ 
quette operators. Besides the On-independent contribu¬ 
tion qq = + |caiP + the conductance contains a 

number of contributions describing quantum interferenee 
of direct tunneling and loop paths {qa^b) or two-loop 
interference {qab)- In a manner detailed below, these 
terms can be used to extract information on the stabi¬ 
lizers Oa,b by measuring the single transport coefficient 
Gi2- 

Surfaee eode manipulation .—Once the code has been 
projected by successive stabilizer measurements to an 
eigenstate {On}, we can manipulate it through the con¬ 
trolled flipping of select plaquettes. To this end, we 
adapt a setup originally proposed for a single Majorana 
wire [ 35 ] and assume that individual MBSs 7 of the code 
are also tunnel-coupled to SETs, i.e., electronic islands 
sufficiently small that the occupancy of a single level can 
be controlled by a nearby gate, see Fig. Describing 
the SET-code coupling by a Hamiltonian as in Eq. 0 
but with Xj X and (0) ^ d, where A is the coupling 
amplitude and d the charge annihilation operator on the 
SET, we consider a situation where the charge occupancy 
of the SET is slowly lowered from 1 to 0, while that of 
a SET coupled to another MBS 7' is raised from 0 to 1. 
This operation amounts to an adiabatic charge pumping 
process during which a single electron enters the code 
through 7 and exits through 7'. (The adiabaticity re¬ 
lies on the pumping rate being slower than the inverse of 
an effective inter-SET coupling amplitude specified be¬ 
low.) The intra-code charge transfer is described by a 
linear superposition of ‘string operators’ 5, i.e., prod¬ 
ucts of tunneling amplitudes 7 zt*^, 7 y establishing a path 
connection between the boxes hosting 7 and 7', respec¬ 
tively. For example, in the situation shown in Fig. the 
shortest string operator is given by S = 7i(2t^2/E^)72- 
Each individual tunneling along the path amounts to a 
virtual excitation of the system, such that a string oper¬ 
ator involving n steps scales ^ E/E'^^ — the process 
is dominated by the shortest path connecting the boxes. 
However, more important in the present context is that 
all string operators commute with the stabilizer system, 
[S, On] = 0, which in turn implies that the action of the 
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Figure 3. Concrete design of a 2D on-chip implementation. 
The Majorana-Cooper box hosts four MBSs (red dots) and 
is implemented as 90 deg rotated H-structure from a pair 
of second-generation topological nanowires m joined by a 
mesoscopic superconducting slab (‘SC’). This slab ensures a 
common charging energy for the box. One normal lead (‘|’) 
and one SET (along with its gate electrode, ‘g’) are tunnel- 
coupled to each MBS by employing T-junctions (see text), 
such that there are eight contacts per box. MBSs on neigh¬ 
boring boxes are connected by short wires via tunnel contacts 

Ct’). 


full tunneling operator, X = A*A' dy Sj y'd'’*', on the 
code eigenstate is effectively described by the pair oper¬ 
ator 77 '. Now each MBS is an element of two diagonally 
adjacent plaquettes, which implies that, depending on 
the separation between 7 and 7 ', a pair ( 7 , 7 ') has either 
two, one, or no plaquettes in common [36]. In the case 
of one shared plaquette, exemplified in Fig. the SET 
operation will create a minimal two stabilizer excitation 
(01,2 in Fig.[2|. For larger separation between the exci¬ 
tation centers 7 , 7 ', four stabilizers get flipped. (Increas¬ 
ing the separation also diminishes the effective inter-SET 
coupling described by X, and this imposes tighter con¬ 
ditions on the adiabaticity of the process m-) Finally, 
successive SET pumping operations can be employed to 
move excitations. 

Operation of the system. — The first step in the oper¬ 
ation of the system containing N plaquettes must be a 
one-time calibration operation in which the non-universal 
but fixed transport coefficients go,AlB,AB appearing in 
the conductance Eq. are determined for all links of 
the system. This step involves the measurement of the 
corresponding tunnel conductances followed by the flip¬ 
ping of one of the neighboring plaquettes [33]. Once these 
coefficients are known, the projective [38] measurement 
of 0{N) suitably chosen tunnel conductances provides 
the full information on the stabilizer system {On}- Ac¬ 
tive operations can then be performed by SET-induced 
changes of stabilizers, by stopping the measurement of 
certain stabilizers, or by the physical removal of bonds, 
cf. Ref. [5]. 

Towards an operational hardware. —Finally, let us sug¬ 
gest a hardware blueprint closer to current-date experi¬ 
mental reality than the schematic drawings in Figs.l^and 

The key ingredient of the setup, single MBSs tunnel- 
coupled to neighboring MBSs and at the same time to 
a normal lead and a SET, could be implemented by so- 
called Majorana T-junctions [39], cf. the junction geome¬ 


try sketched in Fig.[^ Such T-junctions are instrumental 
to Majorana braiding protocols [HHni and a concrete 
roadmap for their realization via crossed nanowires ex¬ 
ists [39]. With the aid of T-junctions, the contact to an 
external probe and the SET associated to the MBS would 
then be implemented nearby and within the 2D structure. 
A four-fold replication of this unit allows contact to the 
four MBSs on a Majorana-Cooper box, realized by the 
small capacitive H-structure in Fig.[^ In a next step, one 
might isolate the entire array of tunnel-coupled boxes via 
a 2D oxide layer, and connect the external (lead and SET 
gate) contacts of each box to electrodes approaching from 
above. The implement ability of a minimal qubit based 
on this idea could be tested on a two-loop prototype with 
six connected boxes, cf. the layout in Fig.|^ which allows 
one to change the system state by SET operations and to 
measure the ensuing change in the tunnel conductance. 
(With minimal changes in the design, one might reduce 
the number of connected boxes to three.) Finally, the 
so-called quasiparticle poisoning time limiting the coher¬ 
ence of MBSs is found to be of order 10 ms [19]. This 
timescale may set a present-date limit to the duration 
of a full operation cycle, i.e., to the measurement of all 
stabilizers. 

Coneluding remarks. —In this Letter, we have sug¬ 
gested a hardware layout and concrete protocols for 
the stabilization and manipulation of Majorana sur¬ 
face codes. Salient features of the proposal include 
implement ability entirely in terms of device technology 
that exists already or is subject to current experimen¬ 
tal effort, a single-step measurement protocol based on 
tunneling conductance measurements, and independence 
of microwave radiation and/or magnetic interferometry. 
Given the current pace of progress, we are optimistic that 
a prototypical qubit based on these ideas might be real¬ 
izable soon. From there, it will still be a long shot to a 
large scale device. However, the absence of external elec¬ 
tromagnetic fields in the operation of the system implies 
a level of protective isolation that can ultimately provide 
a powerful aid towards scalability. 
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In this supplementary material we include a detailed derivation of the effective tunnelling Hamiltonian and 

conductance for the noninvasive measurement. 


We here sketch the derivation of the effective Hamilto¬ 
nian [Eq. (3) in the Letter] and the resulting conductance 
[Eq. (4)1 using the perturbative Schrieffer-Wolff (SW) 
procedure (see Ref. [24] of the Letter for an introduction). 
Eor brevity, let us consider the example of tunneling via 
the ‘direct link’ 7172 between leads 1 and 2, see Eig. 2 
in the Letter. The relevant parts of the Hamiltonian are 
denoted as 

Ht = +h.c. , 

i=i .2 

where jVj and ipj are the charge number operators on 
the two boxes and their conjugate phases, respectively. 
Now consider backgates which are set close to integer 
charges Qj, i.e. Ugj = Qj-\-Angj with \Angj \ 1. The 

energy needed to add (remove) a single charge to the jth 
island is then given by = (1 T We may 

now restrict ourselves to sequences of tunneling which 
implement an effective transfer from lead 2 to lead 1 , i.e., 
^ 4^1 4^2 5 since the other direction follows by Hermitian 
conjugation. The minimal process has to involve each 
of the operators Li ^ ^ 724 ^ 26 *"^^/^, 

and Ai 2 ^ 1"^ once. As an example, the 

sequence L\A\ 2 E\ first brings us to an excited state of 
energy by adding an electron to the box hosting 72 ; 
subsequently, the electron is transferred to the box of 71 , 
yielding the scale E^ ^ last it is removed via lead 1. The 
relevant energy denominator hence follows as 1 /(E]^ ^t)' 
Taking all six possible permutations, corresponding to 
different time ordering of tunneling events, and summing 
up the excitation energy denominators, it is convenient 
to define the ‘detuning parameter’ 77 used in the Letter, 

{Et - E^){Et - E^) 

E+E^E+E^ 

_ 4 4Ang7Ang,2 _ 4 

" Ec,iEc,2 (1 - 4An2 i)(l - 4An2 2 ) "• Ec,iEcA 

The effective tunneling amplitude mediated by the ‘direct 


r 


link’ ^ 7 i 72 , with the tunneling amplitudes in Eq. (S.l), 
then follows as 


Hl2,eS = +h.C., T 


2Ai A2ti2 
Ec,iEc,2^ 


(S.3) 


We observe that in order to obtain a finite tunneling am¬ 
plitude r, we have to detune both backgates n^ 7 =i ^2 away 
from integer values, while Ec,i 7 ^ Ec ,2 has no substan¬ 
tial effect. (We thus assume isotropic Ec below and in 
the main text.) Eollowing the same steps for tunneling 
around the plaquettes Oa,b containing the pair 7172 , 
cf. Eig. 2, we obtain the full effective Hamiltonian as 
stated in Eq. (3) of the Letter. Remarkably, the param¬ 
eter ^ = r/a used in that equation contains no complex 
phase and can be changed through the detuning parame¬ 
ter T]. Note that we have neglected the detuning effect on 
the plaquette contributions, i.e., on the terms ^ Cn ap¬ 
pearing in Eqs. (1) and (3) of the Letter. While this can 
be easily included, it has no significant impact on the 
physics nor on the operation of the code. Calculating 
the tunnel conductance by perturbation theory in 
we obtain Eq. (4) in the Letter. 

A measurement of the conductance G 12 hence contains 
terms proportional to Ob (with coefficient gA,B = 
2 ^Re(c^^B)) and OaOb (with coeff. qab = 2 Re(c^CB)). 
Using a sequence of conductance measurements G 12 and 
controlled plaquette flips on the A or B plaquette, each 
of these values may be extracted separately: the sum of 
two measurements Gi 2 {Ob) and Gi 2 {—Ob) reveals Oa- 
Additionally switching between Oa and —Oa determines 
the corresponding conductance coefficient qa- 
Eor good visibility of conductance switches, we may set 
1^1 ~ Wa,b\ by adjusting the detuning 77 . Assuming we 
are deep in the charge quantization regime, \t\/Ec ^ 1, 
this results in the condition Ang^iAng ^2 ~ 

Eor some \t\/Ec ^1/3 yielding reasonable plaquette en¬ 
ergies, the necessary detuning away from the Coulomb 
valley center can hence be estimated as \Angj\ < 0.16, 
still sufficiently far from half-integer. 








